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Abstract 

We develop the new method of optimized discretization for imaging the rela- 
tive source from two particle correlation functions. In this method, the source 
resolution depends on the relative particle separation and is adjusted to avail- 
able data and their errors. We test the method by restoring assumed pp 
sources and then apply the method to pp and IMF data. In reactions below 
100 MeV/nucleon, significant portions of the sources extend to large distances 
(r > 20 fm). The results from the imaging show the inadequacy of common 
Gaussian source-parametrizations. We establish a simple relation between the 
height of the pp correlation function and the source value at short distances, 
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and between the height and the proton freeze-out phase-space density. 
PACS numbers: 25.70.Z, 25.75.Gz, 25.70.Pq, 29.85,+c 
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I. INTRODUCTION 



Imaging techniques are used in many diverse areas such as geophysics, astronomy, medical 
diagnostics, and police work. The goal of imaging varies widely from determining the density 
of the Earth's interior to reading license plates from blurred photographs to issue speeding 
fines. A typical linear imaging problem involves extracting a source from measured data, 
where the data is represented as the convolution of a kernel with a source. An example of 
this is imaging the relative source function £p(r) from the two-particle correlation function, 
Cp(q), often measured in nuclear reactions [|TJ: 

Cp(q) - 1 = ftp(q) = JdrK(q,r)S P (r). (1) 

Here P is the total momentum of the particle pair and the kernel K is is equal to J2|-Q 

K(q,r) = |$(-)(r)| 2 -l. (2) 

The wavefunction describes the propagation of the pair from a cm separation r out 
to the detector at infinity, where relative momentum q is reached. In the semiclassical 
limit, the relative source Sp for any two particles represents the probability distribution for 
emitting the the two particles with a separation r, in their center of mass. 

On the face of it, imaging appears straightforward: in the case of the two-particle correla- 
tion function, one could just discretize Cp and Sp and invert the resulting matrix equation. 
However, in practice this may not work as small variations in data, within statistical or 
systematic errors, can generate huge changes in the imaged source. This stability problem 
is well known in other fields and vast literature exists on its resolution [f)|-f|. In fact, many 
imaging problems [10, 7j that are now routinely solved would be considered ill-posed in the 



sense of Hadamard, who discussed stability in inversion as early as in 1923 [p]| . A key 



element in stabilizing the inversion process was found by Tikhonov [12] who showed that 
inversion problems can be regularized by using known constraints satisfied by the source. 
Thus, in optical imaging (e.g. in astronomy) the most common method is the maximum 
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entropy method [ I3|, 14 ,[7| ■ The maximum entropy method works by using the most likely 



(i.e. uniform) image at scales below the resolution limit to stabilize the restored image at 
larger scales. In this method, the image is constrained to be positive. In this paper, we 
describe a method for imaging the relative source function from the two-particle correlation 
function, without explicit constraints. Our new method is similar in spirit to the maximum 
entropy method in that we optimize the scale of the kernel K and then use the optimized 
scale to restore the actual source (i.e. the entire image). 

In a recent paper |y, we demonstrated that it is possible to image the relative sources 
in heavy- ion reactions. This represented a step forward as, in the past, the integral relation 
between the correlation functions and the sources was only used to determine the gross 
features of the sources, such as representative radii. In || we put forward two methods 
of imaging sources from correlations. In the first of these methods, we used the fact that, 
when mutual interactions could be ignored or easily corrected such as with photons or pions, 
the kernel is dominated by Bose-Einstein interference. The kernel then becomes a cosine and 
imaging turns to inverting a Fourier cosine transform. The method allowed us to determine 
source values for any separation within a certain range. In the second method, we simply 
discretized the source at equal intervals. This method is more general and can be applied to 
such proton-proton (pp) correlations. We did not discuss the stability of either of 

these methods. 

The question that we now ask is whether discretizing the source at equal intervals is 
optimal in our second, more general, imaging method |l|]. For example, in the pp case, the 
correlation function is dominated by the Coulomb interaction at low-relative momenta and 
by the strong interaction and antisymmetrization at intermediate momenta. The different 
momentum regions should give access to large and short distances within the source, respec- 
tively, with the resolution decreasing at the large distances, rather than being fixed. This 
suggests that, by varying the size of the discretization interval, we can optimize the kernel 
to best restore a given source. Since we do not know the source ahead of time, we must 
specify an "model source" in order to choose the best kernel. This kernel is then used in the 



actual inversion process. 

We illustrate the seriousness of the stability problems by constructing a test correlation 
function, with errors, out of an assumed source and attempting to restore the source. Next, 
we demonstrate how our method stabilizes the inversion. Following this, we apply the new 
method to the pp data JT^,|nj already analyzed in [[!]]. We assess whether the sources are 



positive definite, which would allow for their semiclassical interpretation. We compare the re- 
sults to results from the transport model of reactions |17| . The pp sources exhibit significant 
variation with total pair momentum gate. We investigate the integrals of the source func- 
tions over imaged regions and we show how delicate the analysis of NN correlations can be 
in the different gates. Finally, we apply our method to analyze intermediate-mass-fragment 
(IMF) correlations Hl8| . The IMF sources exhibit significant variation with bombarding 



energy. 



II. IMAGING 

In this section we discuss the general imaging technique. We begin with equations ([[]) 
and (H). Since our goal is to restore the source function Sp, we should discuss briefly 
what it is. The source Sp is proportional to a Wigner transform of a quantal expectation 
value and in the semiclassical limit, the source represents the distribution of relative 



cm separation of the last collisions of the two particles. In this limit, it may be expressed 
in terms of production and absorption rates ||,[[|] and, by nature, it is positive definite. 
Generally, the source is normalized to 1. When pair emission is significantly distorted in the 
vicinity of the source, e.g. due to the Coulomb interaction with the source, those distortions 
may be accounted for either in the wavefunction in the kernel or absorbed into the definition 
of the source. If the latter is the case, or distortions are simply negligible, the functions in ([]]) 
may be expanded into spherical harmonics and relations between angular coefficients in the 
expansion follow [1]]. In particular, the angle-averaged correlation function is related to the 
angle-averaged source through the angle-averaged kernel: 
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C P (g)-l = ft P (g) = 47r / dr r 2 K (q, r) S P (r) . (3) 



o 



The form of the kernel depends on the particle pair. In the pion-pion case the kernel 
is K(q, r) = cos (2q • r), given that the mutual Coulomb interaction is corrected for in Cp. 
The angle-averaged kernel follows simply as Ko(q,r) = sin (2qr)/(2qr). For these kernels, 
the source Sp is an inverse Fourier transform of C [OJ. For the angle-averaged source one 
finds: 

rS P (r) = — / dqq sin (2qr)(C P (q)-l). (4) 
n z Jo 

In practice, the integration in (|]) is cut off at some suitably chosen value of q mauX . This limits, 
from above, the region where the correlation in (||) is dominated by two-pion interference. 
In the proton-proton (pp) case, the angle and spin averaged kernel is U 

K (q,r) = W( 2 3 + 1 ) (flffto) 2 - 1 , (5) 

Z jsW 

where gff a is the radial wave function with outgoing asymptotic angular momentum i. 
In the classical limit of correlations of purely Coulomb origin, such as investigated between 



IMFs pg|, the kernel is 

K {q, r) = 6{r - r c ) (1 - r c /r) 1/2 - 1 , (6) 

with the distance of closest approach r c (q) = 2\iZ\ Z 2 e 2 /q 2 . 

Typically in an experiment, the correlation function Cp is determined at discrete val- 
ues of the magnitude of relative momentum {<7i}j=i,...,M for directionally-averaged function, 
or on a mesh in the momentum space {qi}j=i,...,M when no averaging is done. With each 
determined value C 2 exp some error Oi is associated. It is this set of values {C^ p }i = i t ,„ t M, 
that we use in determining the source function. In the 3-dimensional case we may introduce 
a rectangular mesh in the space of relative particle separation and assume that the source 
function is approximately constant within different cells of the mesh. In the angular expan- 
sion of [|IJ, we may also assume that the spherical expansion coefficients vary slowly within 
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the r-intervals. In the angle-averaged case, discretization amounts simply to the representa- 
tion S ~ J2f=i Sj 9j{ r )i where N is the number of intervals in r, gj(r) = 1 for rj_i < r < Vj, 
and gj(r) = 0, with rj = j Ar. On inserting a discretized form of S into Eq. (|I|) or (Q), we 



find a set of equations for the correlation functions {Ci}i=i,...,M, m terms of {Sj} 



j=l,...,N, 



N 



Cf -l=72f = X;^«5 i> (7) 
where, in the angle-averaged case, 

Kij = 4vr P drr 2 K{q u r). (8) 
The source values may be searched for by a x 2 minimization 

} =nun. (9) 

»=i CT « 

If we do not constrain the space within which we search for {Sj}j=i t ... t N, then a set of 
linear equations for the values follows by differentiating (§) with respect to {S^j^,,,^, 

sr^( K .. Sj -Kr p )K ik = o, (io) 

ij °i 

or in a matrix form 

K T B{KS-TZ cxp ) =0, (11) 
where 5^ = 5u/af. This matrix equation can be solved for S: 

s = (k t bkY 1 k t Bn cxp . (12) 



III. IMAGE STABILITY 

Our discussion of the stability of the imaged source function begins with the errors of 
{5j}j=i,...,jv that we determine from fll2|) , 



^; = (if T M)^ . (13) 



The N x N matrix K B K in (|P2|) is symmetric, positive definite and may be diagonalized, 



M i N 

(K T BK) = E - K ik Kij = E A « < u i > ( 14 ) 

J i=l ^"i o=l 

where {u a } Q =i,...,Ar are orthonormal and X a > 0. With flT3| ) and (0), the square errors for 
individual values of S are 

A 2 ^ = E^- (15) 



We see in fll5|) that the errors of the source diverge (or the inversion problem becomes 
unstable) if one or more of the eigenvalues A approaches zero. In particular, this happens 
when K maps an investigated spatial region to zero. A specific case is when one of of the 
particles is neutral so |$| 2 ~ 1, cf. Eq. (|2|). Moreover, instability can arise when one demands 
too high a resolution for a given set of measurements. In such a case, what might happen is 
that K smoothes out variations in S, so we lose this information in the correlation function. 
If we then try to restore S, we find that we cannot restore S uniquely at high resolution. 
However at lower resolution, we might still be able to restore it. Unlike typical numerical 
methods, we need a singular rather than a smooth kernel |7|,?I[] for our inversion problem to 



be tractable. Finally, a A close to zero can be reached by accident for an unfortunate choice 
of {rk}k=i,...,M for a given measurement. 

To illustrate how serious the issue of stability can be, we take a relative pp source of 
a Gaussian form 

5(r)= (2^)^ eXP ("^)' (16) 
and generate a correlation function C at relative momenta q separated by Aq = 2 MeV/c. 
We use the folding (^) with the wavefunctions in the kernel calculated by solving the 
Schrodinger equation with the regularized Reid soft-core potential REID93 [22j| . This sim- 
ulated source is shown in Fig. [I]. We take Rq = 3.5 fm in the source and we add random 
Gaussian-distributed errors to the correlation function from the folding. The rms mag- 
nitude of the error is 0.015, which is representative of the pp data of Ref. [15| analysed 
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in [|| . We then attempt to restore the source by discretizing it with 7 intervals of fixed size, 
Ar = 2 fm, for r = (0 — 14) fm. We use a g-interval similar to the one used in Ref. |1[, i.e. 
10 MeV/c < q < 86 MeV/c. (Note that a naive Fourier-transform argument suggests that 
a number larger than 7 of r-intervals narrower than 2 fm could be used in the inversion.) 
Results from the restoration according to Eqs. ( |T2"D and ( |TB| ) are shown in Fig. ^, together 
with the original source function (|i~6"l) . Clearly the errors for restored source far exceed the 
original source function. In fact, every second value of the restored source is negative. Inci- 
dentally, this is one of our more fortunate simulations, since all of the errors actually fit on 
the logarithmic plot. 

We next illustrate the dramatic stabilizing effect that the constraints have on the imaging, 



as it was first discussed by Tikhonov ||12|| . We carry out the inversion using the same 
correlation function and errors that we used for Fig. (see Fig. |]). We impose the constraints 
that the imaged source is positive definite, i.e. Sj > 0, as expected in the semiclassical 
limit, and that the source is normalized: within the restored region 1 > An J rjv drr 2 S ~ 
AwYsjLi Sj IrLx drr 2 . Following the general strategies |2B| for estimating values with errors 
and under constraints, we carry out our imaging by sampling the values of the correlation of 
function according to the errors {&i}i=i,...,N (equal to 0.015 in our case) and by applying (P"2"|). 
This amounts to the replacements in (|l~2"l): Cj — > Ci + ai^i, i = 1, . . . , N, with £'s drawn from 
the standard normal distribution. We accept only those samplings where the constraints 
are met. With these samplings, we calculate the average source values and the average 
dispersions. The results are shown in Fig. |3| together with the original source. They now 
compare favorably to the original source. 

In the maximum entropy method 0, as we stated, the constraint of positive source 



definiteness is implicit. Another constraint one might use |12| , p!0| , |2T|] is an assumption of 
smoothness of the source, permitting inversion with more points in the image then there 
are in the data. Cutting off the integral in the Fourier-transform method [1[] at g max implies 
the constraining assumption that S varies slowly on the scale of l/(2g max ) « 1.2 fm in [|TJ, 
which is reasonable given the range of strong interactions. Finally we comment that, in 



imaging terms, the common Gaussian parametrization of sources in heavy-ion collisions is a 
very extreme constraint for stabilizing the inversion. 



IV. OPTIMIZED DISCRETIZATION 

While imposing constraints on the source may stabilize the inversion, we have developed 
an imaging method that can yield very satisfactory results even without any constraints. 
Indeed, one may want to see directly whether the data are consistent with positive definite 
sources. The essential step is to vary the size of the r-intervals to minimize relative error of 
the source. 

Thus, the first stage of analysis involves the values of relative momenta {qi}i=i,...,M, where 
correlation function was measured, and errors on these measurements {crj} i=1 jM , but not 
the values themselves {C i } i= i v .. j M- Specifically, we vary the edges of the intervals for source 

discretization, {fj}^ ^, demanding that the sum of errors relative to some "model source" 

is minimized at fixed N and r = 0, 



N 

E 

3=1 



AS, 



Cmod 



(17) 



where the {ASj}j=i t ,,, t N are the square roots of the {A 2 Sj}j=i t ,, n N in Eq. (|15|). We find 
a rather weak sensitivity of the results to fine details of S mod in (|17D , so we just use a simple 
exponential form S mod oc exp (-r/R' ), S™ od = S^fa-i + r^/2), with R' of the order of 
few fm. The exponential form is consistent with a possible tail in the source due to prolonged 
decays. 

Features of the squared wavefunction in (Q) and the binning in q appear to have the 
greatest affect on determining the best set {r J }j = i r .. ) jv • Nevertheless, it is important to use 
relative errors, with some sensible S mod , in fllTf ). If absolute errors are taken, then the r 2 
weight from angle- averaging in favors large r's. The net result is that we learn that the 
source is close to zero at large r to a very high accuracy; we do not need imaging to tell 
us this. Our practical observation is that the sum of relative errors in fll7D , rather than the 
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sum of squares, is preferred for minimization; the sum of squares pushes {rj} i= i^„^ inwards, 
leaving little resolution at high r. 

Once an optimal set of {rj}i = i t ,.,^ is found, Eq. ( |T2"D can be used to determine the source. 
The results of applying our procedure to the simulated pp correlation function of the pre- 
ceeding section, are shown in Fig. [|, for N = 7. The optimal intervals for discretization 
typically increase in size with r. For example, the first interval in Fig. |] is 2 fm wide and the 
sixth is 3.6 fm wide. The figure clearly shows that we can satisfactorily restore the source 
without imposing any constraints. Figure f| shows the results from a similar restoration of 
the source with an exponential tail: 

1 1 / r 2 \ 1 15 r 

(r) ~ 2 (2^)3/2 ex P y~ ^2 J + 2 4^Rf exptr/flO-l ' { ' 

where R® = 3.5 fm and Ri = 6 fm. We show the correlation function corresponding to the 
restored source in Fig. ([!]), both with and without errors of the rms magnitude of 0.015. 
Since the same N and the same 0i}i=i,...,Af are used in the inversion, we find the same 
optimal {rj }i=i,..., at used in Fig. ||. The restored source gives evidence for the tail in the 
source, despite of the fact that the magnitude of the tail is lower by 2 orders of magnitude 
compared to the maximum at r = 0. Comparing figures |] and [|, we see that our method 
can discriminate between the two source shapes. If we impose additional constraints to the 
optimized discretization method, the agreement between the restored and original source 
functions further improves. 

V. SOURCES FROM PP DATA 

Now we apply our new optimized discretization method to analyse the pp correlation 
data [15,161 , from the 14 N + 27 A1 reaction at 75 MeV/nucleon, imaged in the naive fashion 



in Ref. [0 . Since our new method does not require positive definiteness for stabilization, we 
are able to lift this constraint and verify whether the data favor positive definite sources. 
Further, we do need to normalize the sources to 1 within the imaged region. We also compare 
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the pp sources from data to those from the transport model 17], over a large range of relative 



separations and magnitude of the sources. Past experiences in comparing semiclassical 



transport models to single-particle and correlation data have been mixed p5|JT6|j2^ -p6| for 
this particular reaction and others in this energy range. 



The low relative-momentum pp-correlations have been determined in Ref. [ jT5| , pr5|| for 
pairs emitted around #i a b = 25° from the 14 N + 27 A1 reaction at 75 MeV/nucleon, in three 
intervals of the total momentum: 270-390, 450-780, and 840-1230 MeV/c. The highest lab 
momenta interval corresponds to the highest proton momenta in the participant cm for this 
reaction. These momenta are higher than the average for participant protons and directed 
rather forward. Transport calculations fl|JTT| show that the highest momenta bin is mostly 
populated by pairs from the semicentral to peripheral collisions. The intermediate momenta 
interval corresponds to the magnitude of typical momenta of participant nucleons in the 
forward NN cm hemisphere. The transport calculations show that these pairs stem mainly 
from the semicentral collisions. The lowest lab momenta interval has both participant and 
target spectator contributions. In the latter case, the transport calculations show that pairs 
are mostly from the semicentral to central collisions. According to the transport model, 
the average emission times for protons, from the first contact of the nuclei, in the three 
momenta intervals, from the highest to lowest, are ~ 35, ~ 80, and ~ 110 fm/c, respectively. 

The results of analyzing data using the optimized discretization method are presented 
in Table |. While we quote only the values obtained using wavefunctions for the REID93 
potential, the values for the NIJM2 [^2| differ only by a fraction of l/1000 th . Both the 
results with and without the positivity and normalization constraints are listed in the ta- 
ble. The results are similar, given their errors. The results obtained without constraints 
are generally consistent with positive definite source functions. Only in the highest total- 
momentum interval there might be some indication of a negative minimum in the source 
function. One can also see in Table | that the magnitude of S in the low-r region increases 
by a factor of 2 each time one switches from a lower to a higher total-momentum interval. 

In Fig. ^| we compare the constrained results from the data to the distributions of 
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relative separation of last collision points for protons with similar momentum from the 
transport model PUTTfl. Clearly, the semiclassical model can only yield positive-definite 
source functions. Again, we can see the focussing of the experimental distribution at low r 
as the pair momentum increases. The large-r tails in the distribution at different momenta 
cannot be accommodated with the Gaussian parametrizations used to describe the low-r 



behavior of the sources 15.16 



Generally (Fig. ^|), the Boltzmann-equation model (BEM) yields relative emission point 
distributions that are similar to the imaged data, including the dependence on total pair 
momentum. In fact, the maximae around q = 20 Mev/c (such as in Fig. [I] of the present 
paper) are nearly the same height as the data (see Fig. 2 in f24|). Such findings are somewhat 
surprising for the low and intermediate total momentum intervals. While BEM adequately 
describes high-momentum wide-angle single-particle spectra of protons, which correspond to 



the highest total-momentum interval (see Fig. 1 in [p4|), the model overestimates the single- 
particle proton spectra by as much as 1.5-5 in the two lower momentum intervals ||24| . [] 
Looking closer at Fig. we find that the distributions from data are somewhat sharper at 
low-r in the two lower total-momentum intervals than the the distributions from the model. 
In the next section, we reveal a serious discrepancy when we go beyond a point by point 
examination. 

VI. GENERALIZED CHAOTICITY PARAMETER 

An important quantity characterizing the imaged distribution can be the integral of the 
source over a region where the distribution is significant. We introduce the parameter 

X(r N )= [ drS(r). (19) 

J r<r N 



1 For other comparisons of the transport theory to single-particle data from the same or similar 



reactions see [15,27.25]; overall proton multiplicities are typically overestimated by a factor of 2, 



possibly due to excessive stopping within semiclassical transport in this energy range. 
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that generalizes the the chaoticity A used to parametrize high-energy tttt correlations. 
The term 'chaoticity' stems from early models involving classical currents used to describe 
7r production. The standard chaoticity parameter is defined by fitting the tttt correlation 
function to a Gaussian: 



Ott roo , . 

C(q) - 1 = — / dr r sin (2qr)S(r) ~ Aexp (- 2q 2 R 2 ) . (20) 
q Jo 1 



When one uses this parametrization of the correlation function, one assumes the following 
parametrization of the source at low r: 

5(r) "(2^)^ eXP (~2^)- (21) 



Our chaoticity parameter generalizes the one in (pO ) and (pl|) in two ways. First, because 
our chaoticity parameter is defined in terms of the imaged source function, rather than 
a Gaussian fit to the correlation function, it can apply to any particle pair, not just pion 
pairs. Second it is dependent on the cut-off, r^. The importance of our definition of the 
generalized chaoticity parameter lies in the fact that some particles in the reaction, such as 
pions or protons, can stem from long-lived resonances and be emitted far from any other 
particles. Thus, they contribute primarily to S at large r, possibly outside the region that 
may be imaged. If the kernel K is close to zero or averages out to zero for moderate to 
high q and large r, then the tail in S does not contribute to the deviations of C from 1 in 



Eq. (]T|) or (pop. The imaged region is limited to and may result in A < 1. Access to C 
at lower q may give higher and A closer to 1. However, given experimental limitations 
on measuring C at the lowest relative q, we only really stand a chance to study large r with 
charged particles (and the higher the charge the better). The charged particles with low 
relative q should be emitted far apart, otherwise they acquire relative kinetic energy due to 
mutual repulsion. For the analysed pp data ||15|| , the low-g region is either not available or is 
associated with large systematic errors. Imaging can only extend up to ~ 20 fm. Thus, 
one can expect A(rjv) < 1- In many tttt measurements, e.g. [23], the resolution allows one 



to image jl] regions of comparable sizes, typically 10-20 fm. In contrast to the pp and tttt 
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data, the data on IMF (such as ||18|| ) often extend to low values of relative velocity because 
of the use of fragments with large charges. This permits imaging up to relative separations 
as large as 50 fm. For a discussion of this, see the next section. 

In Table [TT], we have tabulated A(r/vr) in each momentum interval for the constrained 
and unconstrained sources of Section |TTT| as well as sources from the BEM. Comparing 
the results of the BEM to the constrained source, we only find agreement in the highest 
momentum interval. It follows that, when compared to the model, significant portions of 
the source are missing from the imaged regions. This discrepancy is especially pronounced 
in the intermediate-momentum region. Nevertheless, it is comforting that the transport 
model describes the features of the relative source for the high momentum protons since 



it properly describes |24j] the high-momentum single-particle spectra. Now, in BEM no 
IMFs are produced and the IMFs may decay over an extended time, contributing to large 
separations in the relative emission function, as they move away from the reaction region. 
Of course these decays produce some final IMFs, contributing to the relative IMF sources at 
distances similar to those for the pp sources. It may be interesting to see whether a significant 
portion of the relative IMF sources can extend beyond ~ 20 fm, as is apparent for the pp 
sources. We check this in the next section of the paper. 

The disagreements between the data and calculations in both the values of A(tat) and 



the single-particle spectra [24|, for the lower momenta, reveal unphysical features of low- 



momentum proton emission in the transport model. The coarse agreement between the 



measured correlation function and the function calculated using the model in Ref. [24 



in 



the lower total-momentum intervals, mentioned in the last section, is coincidental. Since 
the images show X(tn) < 1, some of the strength of S is shifted out to large r lowering the 
height of C. For BEM, the height of C is lowered because S is softer than S from the data. 
Thus, the BEM correlation function can match the height of the measured correlation peak, 
while not matching the overall shape of the correlation function. For other systems in the 
general energy range, disagreements were found even in the height of C [jn^E^,^] . 

These conclusions raise the question of the general sensibility of attempting to fit the 
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magnitude of pp or nn correlation function at the maximum [^,|T^,|T^] by adjusting the radius 



fitting the tttt correlation functions at intermediate relative momenta, both the strength and 
extent of the source function at low r are varied: A is read off from the magnitude of the 
correlation function at low q and Rq is read off from the width of the correlation function in q. 
However, because of the resonant nature of the low-momentum NN interaction, the height 
of the NN correlation function is determined by the magnitude of the source within the 
resonance peak of the wavefunction. That magnitude is both affected by the overall low-r 
strength of the source and the low-r source falloff. This is illustrated in Fig. [7] which shows 
pp correlation functions for the source (|2T|): the same maximum height can be obtained 
using R = 4.5 fm and A = 1 as using R Q = 3.5 fm and A = 0.5. 

Given that the resonance peak in the 1 Sq wavefunction is quite narrow (it has an outer 
radius of ~ 2.5 fm) and the source falloff cuts off large-r contributions to the integration 
in (fl|), the low-r limit of S is proportional to the C — 1 at the maximum (to a ±20% level) 
for virtually all low-r falloffs that may be encountered in practice, Rq = 2.5 — 6.0 fm: 



For the two sources we used in Fig. [7|, we get about the same value of S(r — > 0) and therefore 
about the same maximum height in C. With the same maximum height, the source falloff 
is reflected in the width of the maximum in Fig. [F[ Note that the value of S at short 
relative distances determines the average freeze-out phase-space density and the entropy in 
the reactions [0,11; the simple relation for estimating the proton phase-space density takes 
the form, from (|22|). 




Rq in the Gaussian parametrization of the source function (see Eq. (|2l|)) with A = 1. When 



C P (20 MeV/c) ^ 1 + 540 fm 3 S P (r 



0). 



(22) 





(23) 



where dN p / dp is the proton momentum distribution. 
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VII. IMF SOURCES 



We now turn to the analysis of IMF sources. We choose the correlation data of Hamilton 
et al. [18], from central 84 Kr + 197 Au reactions at 35, 55, and 70 MeV/nucleon, because 



these data give us the opportunity to examine the variation of sources with beam energy. 

Pairs were collected in the angular range of 25° < 6^ < 50° in order to limit contributions 
to the correlation functions in [I3| from targetlike residues. Hamilton et al. determine the 



functions in terms of the reduced velocity 

v 

Vrcd = (z 1 + z 2 y/2 ' (24) 

under the assumptions that the Coulomb correlation dominated the fragment correlation 
and that the fragments were approximately symmetric, Z/A ~ 1/2. The distance of closest 
approach in (^j) for symmetric fragments is 

Ai A 2 m N v l m N v' ed 

The correlation functions at the three beam energies are shown in Fig. |8|. We have 
reduced the normalization for the correlation function at 35 MeV/nucleon by 5%, compared 



to fll8j[ , in order to better satisfy the condition that C — > 1 at large t> re d, outside of the 
region shown in Fig. In examining Eqs. ([TJ), (||), (ffip, and Fig. ^, certain issues become 
apparent. Contributions to the Coulomb correlation function C for a given t> re d stem only 
from the region of the source with r > r c . As t> re d increases from 0, the distance of closest 
approach r c decreases, with more and more inner regions of the source S contributing to C. 
The low-f re d correlation functions at the three beam energies in Fig. |8| are quite similar, 
suggesting that the tails of the source functions are similar. Differences occur at higher f rec j, 
indicating differences in the inner source regions. 

To image of the IMF sources, we optimize {rj}^^...^ as in the pp case, but we add the 
constraint r\ > r™ in . We do this because the Coulomb interaction in (P does not dominate 
when the measured fragments are in close contact. The Coulomb correlation alone cannot be 
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relied upon to get information on the most inner portion of the source. The typical touching 



distance for the fragments measured in [18|] is r t ~ 5 fm; we chose a value r™ in = 7.0 fm 
which ensures that there is more volume in the lowest bin outside r t , than inside rt- 

The results from the imaging are shown in Fig. ^. Given the errors in the figure, the tails 
of the sources at the three energies are not very different. However, we observe significant 
variation with energy at short relative distances (r < 12 fm) with the source undergoing 
a larger change between 55 and 75 MeV/nucleon, than between 35 and 55 MeV/nucleon. 



In Table |TJ we tabulate the generalized chaoticity parameter for the IMF sources with 
the cutoff of 20 fm and over the whole restored image. The values of A for the whole 
image (cutoff of 90 fm) are all consistent with 1, within errors. The values of A(20 fm) 
20-30% below 1 and Fig. ^| indicate that a large part of IMF emission occurs at distances 
that are not imaged with protons. Nevertheless, the values of A(~20 fm) the low- and 
high-momentum pp sources in Section [V] are roughly equal to the IMF A (20 fm) but the 
intermediate-momentum pp A (20 fm) is lower. It should be mentioned that no complete 
quantitative agreement should be expected, even if the data were from the same reaction 
and pertained to the same particle-velocity range. This is because more protons than IMFs 
can stem from secondary decays. Besides, the velocity gained by a proton in a typical 
decay is large compared to the relevant relative velocities in pp correlations, but the velocity 
gained by an IMF can be quite small on the scale of velocities relevant for the IMF Coulomb 
correlations. Thus, the IMF correlations may reflect the primary parent sources, which are 
concentrated around the origin, rather than the final sources. 

The tails of the IMF sources extend so far that they must be associated with the time 
extension of emission. Even so, it is interesting to ask how far into the center of the source 
must we go to see the effects due to the spatial extent of the primary source. In fI8|| , the 
combinations of single-particle source radii and lifetimes that gave acceptable descriptions 
of their data have radii varying between 5 and 12 fm. In general, the combination of spatial 
extent and lifetime effects should give rise to a bone-like shape of the relative source, with the 
source elongation being due to the emission lifetime. With this, one could try to separate 



the temporal and spatial effects by restoring the full a three-dimensional source. In the 
angle-averaged source, the part dominated by lifetime effects should fall off as an exponential 
divided by the square of the separation, r 2 , as a function of relative separation. On the other 
hand, the part of the relative source dominated by spatial effects may fall off at a slower 
pace or even be constant. In Fig. |S], we see that the sources change weakly with r at 35 and 
55 MeV/nucleon and faster at 70 MeV/nucleon, within the range where sources vary with 
energy (r < 12 fm). For reference, in the insert to Fig. we show the IMF source multiplied 
by r 2 . We see an edge at r ~ 11 fm at 35 and at 55 MeV/nucleon which disappears at 
70 MeV/nucleon, suggesting that the IMF's are emitted from a spatial region, with a radius 
R ~ 11/V2 ~ 8 fm, that becomes more dilute as the energy increases. This observation is 
consistant with the expectation that dilute sources are more effective at emitting IMF's since 
IMF yields in central symmetric collisions or from central sources in asymmetric collisions 
maximize towards 100 MeV/nucleon fl30| . 



VIII. CONCLUSIONS 

We have introduced the new method of optimized discretization to image sources. 
This method is suitable for determining the relative source from different particle pair cor- 
relations in heavy-ion collisions. Recognizing the need for stable imaging, in our method 
we adjust the resolution to minimize the relative errors of the source. The fact that we can 
actually estimate errors in this method gives our method a significant advantage over the 
maximum entropy method. We tested our method by restoring assumed compact pp sources 
and found that the quality of the restored source is comparable to the restored source we 
obtain by imposing the constraints of positivity and normalization. Imposing these con- 
straints in our new method further reduces the source errors, but imposing these constraints 
is no longer required. The new method allows one to study the long-range source structure 
by adjusting the overall size of the imaged region and the resolution at large distances. 

The robustness of our new imaging method gave us the capability to examine the positive 
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definiteness of proton Wigner sources in the 14 N + 27 Al reaction at 75 MeV/nucleon p],[l| 



In principle, unraveling the quantal negative values of a Wigner function is not far fetched. In 



fact, negative values of Wigner functions have been observed in interfering atomic beams 31 



Admittedly, if we had discovered such values in the heavy-ion reactions, we would first have 
to check for a possible breakdown of the assumptions leading to (|I|) and on systematic errors 
in data, before concluding on a success. The extensive averaging in the reactions^ makes it 
unlikely that a genuinely quantal oscillation in the source function would survive except in 
the very tail of the function. 

We found that our imaged sources change significantly with the total pair momentum, 
becoming sharpest for the largest momenta in the cm. Significant portions of the imaged 
source are missing from the imaged regionP] at typical participant momenta in the cm, but 
not at the highest momenta. The chaoticity parameter (the integral of the source) from 
the Boltzmann-equation model |2^JT^] agrees with the data at the highest momenta, but 
the integral is close to one in the participant and target-emission momenta. Nevertheless, 
the model yields the correct height of the maximae of the correlation functions EM . This 
is because the right combination of source normalization and sharpness in the model can 
yield the right value of S at short separations, Sp(r — > 0), and this primarily determines the 



height of C. Gaussian-source fits to the height of the pp correlation function [^ . |i5[|16 



arc 



of a limited value because considerable source strength may lie at large relative separations. 

In our analysis of midrapidity IMF sources in central 84 Kr + 197 Au reactions at dif- 
ferent beam energies, we found a significant variation of the sources with energy at short 
distances, but not at large distances. Considerable portions of the IMF sources extend to 
large distances (r > 20 fm) just like the lower total-momentum pp sources. It would be very 



2 The averaging is over the impact parameter, the central position for the source and emission 
times (R, and t\ and t2, respectively, in Eq. (2) in |fj]), and over the total pair momentum P. 
3 The imaged region corresponds to relative distances with r < 21 fm. 
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interesting to image both the IMF and pp sources in one reaction. There is a deficiency of 
our fragmentation analysis method, namely our method lacks three-body Coulomb effects. 
When weak, these effects could be included as a first-order perturbation. 

As this work was nearing completion, it was suggested |33|] that we use nested Gaussians, 
of variable width and centered at r = 0, in place of sharp-edged intervals in our source 
description. The Gaussians would generalize nicely to three dimensions while maintaining 
the features of the error optimization in our method and we will investigate the possibility 
of using them in the near future. Our strategy of letting the errors and the kernel choose 
what source they can image is novel not just for the problem of inverting correlations but 
the inversion problem in general 0. 
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TABLES 



TABLE I. Relative pp source values restored from the data of Ref. [15] through the optimized 
discretization method, with (S co ) and without (S nc ) constraints imposed for three total momen- 
tum gates. The centers of the discretization intervals are rfc-1/2 = (r& + r&_i)/2 and the size is 
Ark = r k — r k-i- The number of intervals is either N = 6 (two lower momentum gates) or N = 7 
(the highest gate). 
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P-Range 
[MeV/c] 



ffc-l/2 

[fm] 



10 4 x S% c 
[fm" 3 ] 



10 4 x S c k ° 
[fm- 3 ] 



Ar fe 
[fm] 



270-390 



1 

2 
3 
4 
5 
6 



0.9 
3.3 
5.9 
8.7 
12.1 
16.9 



7.6±2.2 
3.55±0.90 
1.07±0.91 
0.85±0.34 
0.30±0.21 
-0.024±0.047 



7.6±2.2 
3.42±0.76 
1.10±0.69 
0.76±0.30 
0.18±0.11 
0.024±0.019 



1.8 
3.0 
2.2 
3.4 
3.4 
6.2 



450-780 



1 

2 
3 
4 
5 
6 



0.9 
3.2 
5.8 
8.3 
11.3 
15.9 



13.89±0.70 
4.18±0.24 
2.69±0.22 
-0.10±0.13 
0.161±0.075 
0.017±0.018 



13.87±0.69 
4.30±0.21 
2.54±0.17 
0.073±0.063 
0.124±0.057 
0.022±0.013 



1.8 
2.8 
2.4 
2.6 
3.4 
5.8 



840-1230 



1 

2 
3 
4 
5 
6 
7 



0.9 
2.8 
5.3 
7.9 
10.6 
13.9 
18.2 



25.0±4.5 
9.8±2.3 

2.89±0.98 
-1.08±0.88 
-0.20±0.38 

0.41±0.22 
0.009±0.089 



24.3±3.7 
11.8±1.7 
2.09±0.59 
0.17±0.16 

0.094±0.088 
0.24±0.10 

0.036±0.030 



1.8 
2.0 
3.0 
2.2 
3.2 
3.4 
5.2 



25 



TABLE II. Comparison of the integral of the relative pp source function, A(rjv), for the restored 
and BEM sources in three total momentum gates in the 14 N + Al reaction at 75 MeV/nucleon. 
The restored sources use the data of Ref. [15|. The integrals are truncated at the distance rjy. 



P-Range A(rjy) 

[MeV/c] unconstrained constrained BEM [fm] 

270-390 0.69±0.22 0.69±0.15 0.98 20.0 

450-780 0.560±0.065 0.574±0.053 0.91 18.8 

840-1230 0.65±0.37 0.87±0.14 0.88 20.8 
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TABLE III. Comparison of the integrals of the midrapidity IMF source function, A(rjv), in 
central 84 Kr + 197 Au reactions at three beam energies, for different truncation points, r^. The 



restored sources use the data of Ref. [18]. 



Beam Energy 


A(90fm) 


A(20fm) 


[MeV/nucleon] 






35 


0.96±0.07 


0.72±0.04 


55 


0.97±0.06 


0.78±0.03 


70 


0.99±0.05 


0.79±0.03 
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FIGURES 




FIG. 1. The solid line is the correlation function from the source in Eq. ( |lq ) and the dashed 
line is from the source in Eq. fll8|). We obtained the wavefunctions in the kernel in (H) by solving 



the Schrodinger equation with the REID93 potential [22|. The symbols are the values of the 
correlation functions with added random errors of rms magnitude 0.015. The square brackets 
above the horizontal axis indicate the range of q we used to restore source. 
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FIG. 2. The solid histogram is the relative pp source function S restored from the simulated 
correlation function in Fig. |l| indicated there by the open symbols. The dashed line is the original 
source function in (|l^) that we used to generate the correlation function. We employed fixed-size 
intervals of Ar = 2 fm for discretizing the source function and we imposed no constraints on S. 
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FIG. 3. This plot is the same as Fig. || except that the restored source is constrained to be 
positive and its integral is restricted to be less than one. 
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FIG. 4. This plot is the same as Figs. || and ||, except that the source is not constrained and 
is restored with the optimized discretization method. 
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FIG. 5. The solid histogram is the relative pp source restored using the optimized discretization 
method. The correlation function used is indicated by the filled circles in Fig. [l|. The original source 
function is shown with the dashed line. 
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FIG. 6. Relative source for protons emitted from the 14 N + 27 Al reaction at 75 MeV/nucleon, 
in the vicinity of #i a b = 25°, in the three total momentum intervals of 270-390 MeV/c (left panel), 
450-780 MeV/c (center panel), and 840-1230 MeV/c (right panel). Solid lines are the source 
values extracted from the data |ll|] and the dotted lines are the source values obtained in the 
Boltzmann-equation calculation. 
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FIG. 7. The solid line is the two-proton correlation function for Rq = 4.5 fm and A = 1 while 
the dashed line is for Rq = 3.5 fm and A = 0.5. The source is the Gaussian in Eq. (^T|). 
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FIG. 8. Fragment-fragment velocity correlation function in central 84 Kr + 197 Au reactions. 
The symbols show the data of Ref. [|18| and the lines show the imaged source function. The 
35 MeV/nucleon data is represented by the open circles and dotted line, the 55 MeV/nucleon data 
by solid circles and dashed line, and 70 MeV/nucleon by open squares and solid line. 
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FIG. 9. Relative source for IMFs emitted from central 84 Kr + 197 Au reactions from the data 
of Ref. |l8| at 35 (dotted line), 55 (dashed line), and 70 MeV/nucleon (solid line). The insert shows 
the source multiplied by r 2 . In both plots, the full image extends out to 90 fm. 



36 



